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Abstract 

Using event-driven molecular dynamics simulations, we quantify how the self diffusivity of con- 
fined hard-sphere fluids depends on the nature of the confining boundaries. We explore systems 
with featureless confining boundaries that treat particle-boundary collisions in different ways and 
also various types of physically (i.e., geometrically) rough boundaries. We show that, for moder- 
ately dense fluids, the ratio of the self diffusivity of a rough wall system to that of an appropriate 
smooth-wall reference system is a linear function of the reciprocal wall separation, with the slope 
depending on the nature of the roughness. We also discuss some simple practical ways to use this 
information to predict confined hard-sphere fluid behavior in different rough-wall systems. 
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I. INTRODUCTION 



Predicting the dynamic properties of moderate-to-high density bulk fluids from first prin- 
ciples remains an outstanding scientific challenge. This already-difficult problem becomes 
considerably more formidable when the fiuid is confined to small spaces, a situation that is 
commonly encountered across a broad range of technologically important settings.'^ In recent 
years, an important step toward solving this problem for confined fiuids was taken by identi- 
fying physically meaningful static (i.e., thermodynamic) variables (or combinations thereof) 
against which dynamic fiuid properties scaled independently of the degree of confinement.'^ 
This identification suggested that dynamic properties such as the self diffusivity of confined 
fluids could be predicted based on knowledge of their static properties and the correspond- 
ing static-dynamics scaling relationship obtained from, e.g., bulk fluid data. To date, tests 
of this scaling strategy for predicting how conflnement affects dynamics of model systems 
have been carried out by molecular simulation of monatomicP^^*^ and a limited number of 
moleculai'^ fluids. Although these tests suggest that the scaling method can successfully 
predict the dynamic behavior of a variety of fluids conflned to pores with different geome- 
tries, little attention has been paid to the nature of the conflning walls themselves and how 
it affects the relationship between static variables and dynamic properties. In this paper, 
we address this issue by investigating how surface roughness impacts the self diffusivity of 
model conflned fluids. 

Smooth, flat structureless walls represent a mathematically convenient and idealized the- 
oretical construct, and thus, historically, have been commonly used in the study of model 
conflned fluids. However, real solid surfaces do in fact exhibit structure, which can sig- 
niflcantly impact the thermodynamic and dynamic properties of the fluids they contact. 
Wetting and lubrication are two obvious examples of phenomena where the structure and 
shape of the fluid-exposed solid surface have signiflcant implications.^^ This should not 
be surprising because surface roughness, which arises from the structural arrangement of 
"wall" particles, ultimately influences the fluid-wall interaction. Also, it is clear that sur- 
face roughness greatly impacts dynamics. For low density gases, where Knudsen diffusiorP' 
dominates, self diffusion decreases with surface roughness.'^^'^ Furthermore, experiments of 
conflned colloids flnd that some particles stick to walls, further slowing dynamics. Simula- 
tions have shown that position-dependent relaxation processes near rough surfaces are much 
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slower than those near smooth surfaces .'^^'^ 

A question which has recieved comparatively less attention is how surface roughness im- 
pact the self-diffusivity of moderate-to-dense confined fluids. The goal of this paper is to an- 
swer this question systematically via molecular simulations. The most obvious way in which 
roughness influences dynamic properties is through the modification of fluid particle-wall 
collisions. Here, we investigate how different representations of this collisional modification 
lead to changes in the self-diffusion of a simple, confined fluid. As a starting point, we study 
the monodisperse hard-sphere (HS) fluid confined between smooth hard walls. Despite its 
simplicity, we choose the HS fluid because it captures many of the important effects arising 
from the dominant excluded volume interactions in liquid^^^^ Obviously, a smooth hard 
wall does not possess any roughness whatsoever, but it still serves as a useful reference point 
for studying confined fluids, specifically the dynamic properties in this paper. In fact, the slit 
pore geometry is often assumed when determining the pore size distribution from adsorption 
isotherms in real porous materials.'^^'^ In our simulations, we introduce surface roughness 
in two separate ways (i) by modifying the boundary conditions for collisions between fluid 
particles and a flat wall surface and (ii) by giving the surface of the confining walls physi- 
cal roughness or shape. The former crudely represents surface roughness on a length scale 
much smaller than the diameters of the fluid particles (e.g., confined colloids), while the 
latter models physical roughness on a length scale comparable to the diameter of the fluid 
particles (e.g., confined molecular fluids). While neither of the above types of walls can be 
said to completely represent the roughness present in real physical systems, they represent 
the types of roughness that can be incorporated in molecular simulations, and serve as a 
good starting point to address the impact of surface roughness on self-diffusivity. 

We find that surface roughness reduces average particle mobility relative to a smooth flat 
wall, where the collisions between fluid particles and the wall surface are perfectly reflecting. 
Moreover, the reduction appears to be systematic with increasing degree of confinement (i.e., 
decreasing pore width). In the case of physically rough walls, we show that it is necessary 
to make the distinction between the spatially homogeneous and inhomogeneous directions 
because the associated self-diffusion coefficients can differ significantly. Finally, we find that 
the dynamic properties of the hard-sphere fluid confined between rough surfaces can be 
regarded as a perturbation on the dynamic properties of an appropriately chosen smooth- 
wall reference system. 
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This paper is organized as follows. In Section ^ we describe the model fluid studied 
in this work, as well as the simulation methods used. We then discuss how wall-surface 



roughness is implemented in Section |III[ Results are presented and discussed in Section [IV 
Finally, we present conclusions and directions for future work in Section |V| 



II. MODEL FLUID 



Event-driven microcanonical molecular dynamics (MD) simulation^ were used to study 
a fluid composed of hard spheres of diameter a and mass m with interaction potential 



oo rij < a 



r^,)=< (1) 
> a, 

where is the distance between particles i and j. Event-driven MD consists of four basic 
steps: (1) Calculate future events (collisions) times. (2) Sort events to determine next event 
to occur. (3) Advance system to next event using Newtonian dynamics (free flight). (4) 
Execute event (collision) to determine new particle velocities, and return to (1). Various 
methods are available to speed up the general algorithm. The exact method used in this 
study is described in Refj28l A rectangular simulation cell of dimensions Lx 
with = 4000 particles was used throughout this work. Periodic boundary conditions 
were applied in the x and y directions in all cases, and in the z direction for the bulk 
fluid. For the confined fluid, where the z direction was non-periodic, appropriate wall- 
boundary conditions (see below) were employed. The dimensions of the simulation box 
were chosen to correspond to a desired reduced number density pa^ (or packing fraction= 
TTpa^/6), and such that the dimensions of the simulation box in the periodic directions were 
(nearly) equal. Initial configurations for the MD simulations were generated by randomly 
inserting hard spheres at low density with no overlaps, followed by compression and (Monte- 
Carlo) relaxation steps until the desired density was obtained. The systems were deemed 
equilibrated when the number of collisions (particle-particle and particle-wall) per unit time 
reached a constant value as a function of time. Equilibration times of 1000 (ma^/fcBT")^/^ 
were found to be sufficient to meet this requirement. Production simulations spanned times 
exceeding 10^ (ma'^ /k-QTY^'^, which allowed the slowest system to display displacements on 
the order of SOcx. The (pore-averaged) self-diffusion coefficients were calculated by fitting the 
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long-time behavior of the mean-squared displacement in a periodic direction to the Einstein 
relation, e.g., = limt_5.oo (Aa;^)/(2t). Where appropriate, self-diffusion coefficients in 
equivalent periodic directions were averaged together. 

III. SURFACE ROUGHNESS AND BOUNDARY CONDITIONS 

We study the hard-sphere fluid confined between parallel hard walls in a slit-pore ge- 
ometry. To avoid any confusion, we point out that none of the walls studied here exhibits 
internal structure. That is, the walls are of uniform density and can be regarded as a con- 
tinuum solid. Let H denote the average distance between the two surfaces of the confining 
walls. In this work, we study two general types of wall surfaces. The first involves fiat walls 
where the boundary conditions for the particle-wall collisions have been modified. These 
walls model surface roughness on a length scale much smaller than the diameters of fluid 
particles, and thus, because the surfaces are geometrically fiat, they are referred to as fea- 
tureless surfaces (walls). The second type of surface studied here possesses physical (i.e., 
geometric) roughness, i.e., the height of the wall surface varies with lateral position, and 
thus they are referred to as physically rough surfaces (walls). Schematics of each type of 
surface are given in Fig. [TJ 

A. Featureless Walls 

We considered three different types of boundary conditions for particle collisions with fea- 
tureless walls. As a useful reference state, we first considered fiat, smooth walls (SW), where 
the fluid-wall collisions are specular. The second type of featureless surface studied was a 
so-called thermal wall (TW). This model boundary attempts to reproduce particle collisions 
with a surface uneven on length scales considerably smaller than the particle diameter, the 
essential feature being that pre- and post-wall collision velocities are uncorrelated. More 
importantly, the distribution of post- collision velocities is determined by the temperature 
of the wall. The third type of featureless surface gives an alternative approach to generate 
uncorrelated pre- and post-wall-collision velocities. We refer to it as a rotational wall (RW) 
for reasons that will become apparent below. 

In order to distinguish between the featureless walls with different boundary conditions. 



5 



consider a particle with a pre-wall-coUision velocity v = {vx, Vy, Vz), and a post-wall collision 
velocity v' = {v'^,v'y,v'^). For the cases considered here, the components of v' are as follows. 

• Smooth walls (SW): v'^ = Vx, v'y = Vy, and v'^ = —Vz 

• Thermal walls (TW): v'^ is chosen according to the following distribution!^ depending 
on wall temperature T„ (set equal to fluid temperature T): 

where the sign of v'^ depends on the location of the wall collision (i.e., positive for lower 
wall collisions and negative for upper wall collisions), v'^ (where i = x,y) is chosen 
randomly from a Gaussian distribution: 



M^z) = jr^ exp 



(2) 



^M) = \ L Trp exp 



11 ~\ 

—mvf 



2kBT^ 



(3) 



Rotational walls (RW): v'^ = —Vz, and the lateral components v'^ and v'y are determined 
by rotating the corresponding components of the pre-wall-coUision velocity Vx and Vy 
by an angle 6: 



v'^ = Vx COS 9' — Vy sin 9' 



v'y = Vx sin 9' + Vy cos 9' , 



(4) 



where 9' = ±9, with the sign chosen randomly. Note that = corresponds to smooth 
walls, while 9 = tt corresponds to bounce-back boundary conditions. 

Note that the boundary conditions described above affect the dynamics of the system 
only. Because all three of the above walls have the same geometrical form and produce the 
same velocity component distributions, they lead to identical thermodynamic and structural 
propertied^. This was verified by the simulation results. 

We note that the featureless walls described above are not intended to model any specific 
physical system. Rather, they are mathematically convenient course-grained models for 
surface roughness. For example, the motivation for the thermal walls goes back to MaxwelP^, 
who considered collisions with a highly uneven low density granular surface. Particles that 
strike this type of surface undergo a series of collisions with many different surface molecules. 
The resulting outgoing velocity is expected to be randomized, with a distribution determined 
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by the temperature of the wall. We also stress that the rotational walls introduced here 
are not intended to mimic any type of real system. Instead, they should be considered 
a mathematical construct that allows us to change the surface boundary condition in a 
continuous fashion from perfectly smooth (9 = 0) to rough (bounce-back, 6 = n) walls, with 
the system having identical thermodynamics for all values of 6. This allows us to study 
how different levels of surface roughness impact the dynamics of the model fluid without 
changing thermodynamics. 



B. Physically rough walls 

For a physically rough surface, we intuitively expect that the height of the surface should 
vary with lateral position. For simplicity, we use the following expressions for the upper and 
lower wall surfaces in a slit- pore geometry: 



My) 
My) 



H 



+ ttw COS 



'27ry 
X 



H 



a,,, cos 



A 



(5) 



where the subscripts L and U denote the lower and upper surfaces, respectively, and 
and A are the amplitude and wavelength of the well-behaved surface variations, respectively. 



Figure [Tb] displays the geometry of the physically rough walls employed. As written in 
Eq. (jsj) and depicted in Figure lb , minima in the upper wall and maxima in the lower wall 
are aligned, which imposes a maximum value upon if a fluid particle is to have access to 
the entire pore length. The average surface-to-surface distance is H as long as the quantity 
Lyl A is an integer, which is true throughout this work. Because we have chosen to have the 
height of the surface be a function of only y, the system is spatially inhomogeneous in the 
y-direction and homogeneous in the x-direction. This will have important consequences for 
the self diffusivity of a fluid confined between these surfaces. 

The rigid character of the confining walls was maintained by requiring perfectly reflecting 
specular particle-surface collisions. However, instead of solving the set of nonlinear equations 
to determine particle-surface collision times, which is computationally expensive, the curved 
walls were discretized into a set of short, connected line segments of length 0.015(j. Using 
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smaller segments does not lead to noticeable changes in the resulting properties of the 
systems studied here. 

We stress that the specific form of the physically rough walls studied here is not intended 
to mimic a specific physical system. Rather, they allow us to systematically determine the 
impact of different surface feature sizes on the self-diffusivity of our model fluid. 

IV. RESULTS 

A. Featureless Flat Walls 

We flrst present results for the hard-sphere fluid conflned between smooth, flat hard 
walls. In Fig. [2| we plot the self-diffusion coefficient, D^^, where the superscript SW 
signifies smooth wall function of the total density p = N/{AH). Data are shown 

for a number of pore widths H. As previously pointed out elsewherd^, all data points seem 
to fall approximately onto a single curve independently of H for pcr^ < 0.75. Below, we 
investigate how modifying the nature of the confining surfaces changes the results obtained 
in this basic reference system. 

B. Thermal Walls 

Figure [s] displays the self-diffusion coefficient for the hard-sphere fiuid confined 

between featureless thermal walls as a function of p for various values of H. Compared to 
the smooth, fiat wall, roughness due to the thermal wall has a clear and noticeable infiuence 
on particle dynamics. Thermal walls reduce particle mobility, with the magnitude of the 
reduction increasing with decreasing H. That average mobility of fiuid particles confined 
between thermal walls decreases with p at fixed H has been previously reported.'^ 

The interesting trend we quantify in this study is the reduction of D'^^ with increasing 
degree of confinement at fixed p. This behavior can be understood by considering the 
quantity the fraction of collisions in the system involving the surfaces of the confining 
walls. Intuitively, we expect non-specular particle-surface collisions, such as those that take 
place in the thermal-wall systems, to slow dynamics in the transverse direction relative to 
specular particle-surface collisions. We base this expectation by considering the free fiight of 
a single particle. In this case, specular particle- wall collisions do not change the transverse 
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motion, while non-specular particle-wall collisions do. We also expect that the fraction of 
wall collisions grows with the prominence of the walls (the fraction of the fluid particles 
near the walls), that is, with decreasing H. That is, in order to maintain a fixed and 
V , a decrease in H must be compensated by an increase in fluid-exposed wall surface area. 
Figure |4] shows that this is indeed the case. Decreasing H at constant p systematically 
increases the fraction of wall collisions. Clearly, the reduction in D^"^ with decreasing H 
observed in Fig. [3] is directly linked to the increasing fraction of wall collisions. Also, note 
that the results given in Fig. |4] are not unique to the thermal- wall systems. In fact, all three 
kinds of featureless surfaces studied in this work yield identical wall-collision statistics. As 
noted above, the collision boundary conditions only affect the dynamic properties of the 
fluid and not the thermodynamic properties. The thermodynamic pressure of the system is 
intimately related to the wall-collision statistics. 

While there are clear differences between the smooth-wall and thermal-wall self- diffusion 
coefficients, and D'^^ , respectively, these differences appear to depend systematically 
on H. This suggests that it might be possible to develop an approach to predict the self- 
diffusivity of the hard-sphere fluid conflned between thermal walls using a limited amount of 
information. In particular, note the dependence of the fraction of wall collisions on density. 
Initially, increasing density leads to a pronounced decrease in the fraction of wall collisions. 
This is due to the associated increase in particle-particle collisions compared to particle-wall 
collisions (not shown). However, the fraction of wall collisions eventually becomes a weak 
function of density. Since we expect the reduction in self diffusivity due to the presence of 
rough walls to scale with the fraction of wall collisions, we should likewise expect a similar 
density dependence of the ratio D^"^ / . 



In Fig. 5a, we plot the quantity /^tw^^sw ^ function of density. The ratio of the 
diffusivities is taken at the same thermodynamic state of the fluid, namely at the same 
density p and average wall separation H. Notice that, for each pore width, the ratio of 
self-diffusion coefficients initially increases with density before reaching a limiting value. 
Furthermore, the behavior of D^"^ / vs p reflects that of fyj vs p. Again, this points to the 
strong connection between the self-diffusion in a system with rough walls to the prominence 
of the walls (i.e., the fraction of the fluid near the walls). The diffusivity ratio takes its 
smallest value at low density, conditions where surface collisions are most influential, and 
thus, where the greatest difference between the two surfaces is observed. Also, the density 
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at which the ratio reaches a plateau value (pcr^ ^ 0.2) is independent of pore width. This is 
due to the dominant influence of particle-particle collisions on the mobility of hard-spheres 
at moderate-to-high densities. 



The most striking feature of Fig. 5a involves the if-dependent plateau values for pa"^ > 



0.2. In Fig. 5b, we plot all of the data points for which pa^ > 0.2 versus the inverse pore 



width H ^. Fig. 5b shows that the ratio of diffusivities can be well described in this density 
range by a linear function of 

^ C T,), (6) 



where the constant C is a fitting parameter, _D™"s'^ is the self-diffusivity of the fluid between 
the rough surface of interest, and D'^^^ is the self diffusivity of the fluid at an appropriately 
chosen reference state under the same thermodynamic conditions. We note that Eq. ([6]) has, 
to first order, the same dependence on H as the self-diffusivity of a Brownian particle in the 
center of a slit-pore.^^ This, in fact, was the inspiration for the functional form chosen in 
Eq. g. For the thermal walls, we take = jjrd ^ ^sw^ ^^^^ 

C ^ 1.44(T. 

Assuming sufficient smooth-wall self-diffusivity data are available as a function of den- 
sity, this provides the basis for estimating the diffusivity of the hard-sphere fluid confined 
between thermal walls for pa^ > 0.2. To estimate the diffusivity of a fluid confined between 
thermal walls in a slit pore of width H at density p, one calculates two quantities. The first 
quantity is the if-dependent ratio £)TW^£)SW ^gjj^g jg| g^n^j second quantity is the 
self-diffusion coefficient of the reference system, Z)^^, which, as discussed in Section [l| may 
be approximated by a scaling analysis. Knowledge of these two quantities then allows for an 
estimate of the self diffusivity between thermal walls. In Fig. [6} we test the predictive ability 
of this approach. The ratio of the predicted to observed thermal wall self diffusivity is plot- 
ted against the observed self diffusivity for the thermal-wall system. If the predictions were 
perfect, the points would fall on the horizontal line /}TW^£)SW _ ^-^^ ^j^jg jg clearly not 
the case. However, the predicted data points do fall within the indicated 5% error bounds. 
Moreover, we emphasize that to truly make predictions, one must still have knowledge of C 
and the reference self-diffusivity. 
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C. Rotational walls 



In Figure [?[ we present results for the self-diffusion coefficient of the hard-sphere 

fluid confined between rotational walls as a function of density p with (a) 6 = 7i/2 and 
various values of H, and (b) H/a = 7.0 and various values of 6. Qualitatively similar 
results were observed for other choices of rotational-wall parameters and can be found in 
Supplementary Materials.'^Sl Recall that 6 is the angle of rotation that the lateral velocity 



of a particle undergoes after it collides with the surface. Fig. 7a shows that self diffusivity 
decreases with increasing density at fixed H and fixed 6. In addition, at fixed density p, the 
self-diffusivity decreases with decreasing H. This latter trend, which was also observed in the 
the thermal-wall systems, is not surprising considering rotational walls, like thermal walls, 
inherently retard particle mobility, an effect that is most apparent at small wall separations 
(at the same fluid density). Therefore, recalling that the collision statistics are independent 
of the boundary condition for featureless walls, the same physics explaining the reduction 



in mobility due to thermal walls mentioned in Section IV B also applies to rotational walls. 
We also expect that the reduction in mobility should increase with 6, since 6 controls the 



effective roughness of the walls. The data presented Figure 7b bear this out. 



Figures 8a and 8b again show that rotational walls influence self diffusivity relative to 
smooth walls in a manner qualitatively similar to thermal walls. Specifically, the ratio of self 
diffusion coefficients /)RW^£)SW initially increases with density and then levels off beyond 
pa^ > 0.2. In fact, for a given 9, the plateau value appears to be systematically dependent 
upon H, which suggests that an approach similar to that adopted for the thermal wall system 



can be used to estimate the diffusivity in the rotational-wall system. Figure 8b shows that 
the ratio /^RW^^sw (iecreases with increasing 9, indicating that the constant C in Eq. ^ 



is dependent upon 9. This just reflects the 9 dependence noted in Fig. 7b Accounting for 
this 9 dependence. Fig. 8c shows that diffusivity ratios with densities pa^ > 0.2 are well 
described by Eq. ([6]). The inset to Fig. 8c displays the values of C{9) from the fits. For 
simplicity, we fit C to the function C = 2.0sin{9 /2), which as shown in Fig. 8c fits the data 
well, specific form of the fit was based on the fact that C = for = (by definition), C is a 
maximum at ^ = vr, and C is periodic in 9. Using the procedure outlined above for thermal 
walls, the self-diffusion of the hard-sphere fluid between rotational walls can be estimated 
in an analogous fashion. Figure |9] tests these estimates, plotting the ratio -Dp^/D^^ versus 
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D , where D^^^ is the predicted value and D is the value observed from simulation. In 
the case of rotational walls, the estimation procedure yields predictions that fall within 10% 
of the actual values. More accurate predictions can clearly be obtained by making use of a 
more quantitative model for C{6). 



D. Physically rough walls 

Here, we present results for physically rough surfaces at fixed wavelength X/a = 3 and 
several densities p = N/V, where V is the surface accessible volume, while (1) fixing the 
amplitude of the features a^/a = 1.0 and varying H and (2) fixing H/a = 7.5 and varying 
ttw Since the surface height of a physically rough wall depends on lateral position, these 



parameters (see Eq. |5]and Fig. lb) allow us to study the infiuence of surface height variation 
on mobility in a systematic way. Results for A/a = 1.5, 2.0, and 6.0 give results qualitatively 
similar to those presented below, and can be found in Supplementary Materials.'^ 



Notice that in the geometry of the physically rough wall system studied (see Fig. lb), the 
X and y directions, though both infinite, are not equivalent. Because the surface height is a 
function of y only, the y direction is termed rough (inhomogeneous) , while the x direction 
is termed smooth (homogeneous). Also, because of this height variation, we expect the self 
diffusivities in the x and y directions at the same state point to be significantly different. 



with Dx > Dy. Figure 10 shows these two self-diffusion coefficients separately as a function 
of density for various values of H and a^. Observe that the self diffusivity in the x direction 
has little dependence on H or a^. That is, the D^-p correlation is almost equivalent to the 
bulk correlation, much like the case for hard spheres between smooth flat walls (see Fig. |2]). 
This seems consistent with previous results for hard-spheres in cylindrical pores where the 
surface is smooth and exhibits curvature.!^ 

In contrast, the self-diffusion coefficient in the y direction depends strongly on the degree 



of confinement (Fig. 10b) and on the amplitude of the surface variation (Fig. lOd). Specif- 
ically, decreasing H at fixed wall feature size (constant a^} and A) systematically decreases 
Dy. This behavior resembles that for the fiuid confined between thermal walls (Fig. [3]) and 
rotational walls (Figs. 7a). Also, at fixed if, Dy systematically decreases with increasing 



(i.e., increasing surface roughness). This effect is analogous to increasing surface roughness 



{6) in the rotational wall system (see Fig. 7b). Also, the slowing down of dynamics due to 
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physical roughness is consistent with previous studies-l^^ 

The disparity between and Dy grows with decreasing H at fixed and A, and with 
increasing roughness {a.^) and fixed H. Given the close parallels between Dy and the self- 
diffusion coefficient of the ffuid between featureless rough walls, we expect the reduction in 
Dy due to surface roughness should be connected to the fraction of collisions in the system 



involving the walls. In Figure 11, we plot fw versus density for the state points considered 



in Fig. 10, Fig. 11a shows that grows with decreasing H at fixed p and a^, while Fig. lib 



shows that the wall-collision fraction grows with at fixed p and H. In both cases, 
grows with the increasing prominence of the walls (caused by decreasing if, or increasing 
ttw) Also, at fixed wall conditions (e.g., and A), initially decreases with increasing 
density (< 0.2). However, this quantity changes only moderately with further increases in 
density (> 0.2). This behavior is analogous to that observed in the featureless-wall systems. 

Because of the qualitative similarities between the properties of the physically rough-wall 
system and the featureless rough-wall systems studied above, we now examine the diffusivity 
ratio Dy/D^'^^. However, unlike the other surfaces encountered in this work, the choice of 
D"^^^ is not obvious for the physically rough walls. For the thermal and rotational walls, we 
chose the smooth (ffat) wall self-diffusion coefficient at the same p and if, which corresponds 
to the self-diffusion of a ffuid at the same thermodynamic state without surface roughness. 
Likewise, for physically rough-wall system, the reference state should be the self-diffusion 
of a system at the same thermodynamic state, but without surface roughness. For this, 



we choose D'^'^^ = D^. In Figure 12a, we plot the ratio Dy/Dj. versus density for various 
values of H and a^. At fixed ii, the ratio initially increases with p and then levels off 
beyond pa > 0.2. For a specified value of p, the ratio DyjD^ decreases with decreasing H 
at fixed surface roughness (Fig. |12a ) and with increasing roughness at fixed H (Fig. 12b). 
Comparison with Figs. 5a and |8] shows that physical roughness alters Dy relative to D^ in 
a manner similar to how featureless rough walls alter the self-diffusion relative to smooth 



walls. Figure 12c shows that, for given surface features, Dy/D^ is a linear function of 
H-^. That is, Dy follows Eq. (g with £1'^°"^^ = Dy and D'""^ = D^. Figure [l2d| shows 
the fit parameter C versus for the different surface feature wavelengths A studied (see 
Supplementary material^. Clearly, C is much more sensitive to than A. We find that a 
linear fit to C as a function of a^, although crude, describes that data reasonably well. 
From the above analysis, we can now formulate a means to predict the self-diffusion in 
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the rough y direction from the self diffusion in the smooth x direction. Specifically, from 



Eq. ([6]), D^^'^'^ = D,j.[l — C (H/a)^^]. Figure 13 shows the quality of the prediction based on 
this formalism for the state points in Fig. 10 In Figs. 13a[ |c we use C from the fits at a given 
A (i.e., data points in Fig. 12d), while in Figs. 13bp we use the approximation C = 2.12a^/(j 
(linear fit in Fig. 12d). We find that using the C for a given wall configuration yields good 
predictions, with the majority of the data within 10% of the actual data. On the other hand, 
using the approximate C yields predictions within 25% of the actual value. The information 
necessary to estimate Dy in this case is greatly reduced, but one would still need knowledge 
of to truly make predictions of Dy. 



V. CONCLUSIONS 



We have systematically studied how surface roughness affects the self-diffusion of confined 
hard-sphere fluids. Specifically, for pa^ > 0.2, we have shown that the ratio of the self- 
diffusion coefficient in a rough wall system to the self-diffusion coefficient in smooth-wall 
reference system is a linear function of reciprocal wall separation [see Eq. (|6|], with the 
precise slope depending on the specific nature of the surface roughness. If this slope and the 
reference self-diffusivity are known, accurate predictions of the self-diffusion in rough wall 
systems can be obtained. 

In real world applications, however, the specific nature of surface roughness is often un- 
known. This study leads to two possible ways of making useful self-diffusivity estimates 
without this knowledge. First, our analysis shows that a only a limited number of mea- 
surements are needed to determine C. Second, if only qualitative information is sought, 
our analysis shows that C is typically of 0(1), and that the basic physics of self-diffusion 
in a rough system relative to that in a smooth system is captured by the correction 
in Eq. ([6]). That is, the specific nature of the surface roughness does not appear to impact 
radically or significantly the resulting dynamics. This is particularly important since mod- 
eling the specific details of real surfaces is a daunting task. However, it is not our intention 
to suggest that any type of surface roughness can be modeled in any way one chooses, for 
example, modeling geometric roughness with featureless rough walls. Rather, the specific 
details of the surface roughness do not appear to be overly important. 

To make predictions for self-diffusion of confined fluids between rough surfaces using 
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Eq. (|6]), one also requires a prediction of the reference state self-diffusion. As discussed in 
Section [T| one may use a scaling method to map the bulk system to the confined system. For 
the case of smooth featureless walls, for example, this has been shown to be very effective.'^ 
However, for the case of the physically rough system, one must be able to predict D^- We 
plan to address this in a future publication. 

Another question is whether the results presented in this paper apply to more complex 
pore scenarios, such as diffusion through porous materials. Such systems have complex ge- 
ometry and connectivity, as well as rough surfaces. We believe that the simplified geometry 
studied here, slit pores, is applicable, however, when one considers that many of the prop- 
erties in such complex systems can be considered as an average over a collection of simple 
slit pore (or cylindrical pore) systems.'^^'^ 

The lack of particle-particle or particle-wall attractions is an obvious shortcoming of this 
study. Clearly, wall attractions will greatly alter the number of particles near the wall, 
and therefore the fraction of wall collisions. As we point out, the ratio of self diffusivity in 
a rough-wall system to that in a corresponding smooth-wall system appears to be closely 
linked with the fraction of wall collisions. Therefore, attractions can have a strong effect 
on the behavior of the self diffusivity ratio. However, we speculate that the introduction 
of moderately strong wall attractions to the rough systems studied here will result in an 
increased retardation of self-diffusion and not a fundamental change to the physics. We also 
plan to explore this issue in detail in a future publication. 
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Figure Captions 




FIG. 1: Schematics of the geometries of the confined fluid system, as described in the text, (a) 
Featureless walls and (b) physically rough walls. 




FIG. 2: Self diffusivity D versus total density p = N/{AH) for the hard-sphere fluid conflned 
between smooth walls in slit-pore geometry. 




FIG. 3: Self diffusivity D versus total density p for the hard-sphere fluid confined between thermal 
walls, for various wall separations H. 
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FIG. 4: Fraction of collisions in the system involving the wall surface fw for the hard-sphere fluid 
confined to a slit-pore geometry. Symbols have the same meaning as in Fig. |3] 
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FIG. 12: Ratio of self diffusivity in rough direction Dy and smooth direction Dx with A/cr = 3.0 
versus (a,b) p and (c) H^^ for the hard-sphere fiuid confined between physically rough walls. 
Symbols in (a-c) are the same as in Fig. 10 In (c) all data points are for pa^ > 0.2 and lines are 



fits to the form Dy/Dx = 1 — C{H/a) ^. Panel (d) displays the values of C versus a^j at all the 
wavelengths studied. 
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